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rN Abstract 

We study the constraints of crossing symmetry and unitarity in general 3D Confor- 
mal Field Theories. In doing so we derive new results for conformal blocks appearing 
in four-point functions of scalars and present an efficient method for their computation 
in arbitrary space-time dimension. Comparing the resulting bounds on operator 
dimensions and OPE coefficients in 3D to known results, we find that the 3D Ising 
model lies at a corner point on the boundary of the allowed parameter space. We also 
derive general upper bounds on the dimensions of higher spin operators, relevant in 
the context of theories with weakly broken higher spin symmetries. 
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1 Introduction 

This paper is the first in a series of works which will conceivably lead to a solution of the 
Conformal Field Theory (CFT) describing the three dimensional (3D) Ising model at the 
critical temperature. Second-order phase transitions in a number of real-world systems are 
known to belong to the same universality class: most notably liquid-vapor transitions and 
transitions in binary fluids and uniaxial magnets. 

Field-theoretical descriptions of critical phenomena and computations of critical expo- 
nents have a long tradition [1]. One well-known approach to this problem is the e-expansion 
[2]. In this method the critical exponents are computed using the usual, perturbative 
field theory in D = 4 — e dimensions, and the physically interesting case of D = 3 is 
obtained by extrapolating to e = 1. The obtained series in e are divergent and need to be 
resummed. Apart from small ambiguities, the final results for the critical exponents agree 
well with experiments and with a host of other approximation techniques (high-temperature 
expansion, Monte-Carlo simulations, etc.). 

In this paper we will develop an alternative method for determining critical exponents 
in D = 3, based on Polyakov's hypothesis of conformal invariance of critical fluctuations 
[3], which was a major motivation for the development of Conformal Field Theory. CFT 
methods have been extremely fruitful in D = 2, allowing one to solve many models of 
critical behavior [4]. The novelty of our project is to apply them in D = 3. The existing 
quantitative approaches to critical phenomena in D = 3 do not take full advantage of 
conformal invariance. 

The CFT describing the 3D Ising model at criticality is not known to possess any 
additional symmetry apart from conformal invariance and Z2 invariance. For this reason we 
will be able to rely only on the most general properties of conformal theories. The study of 
such general properties goes back to the 1970s. The required fundamental concepts are the 
classification of primary operators, the conformally-invariant operator product expansion, 
conformal blocks, and the idea of the nonperturbative conformal bootstrap, which were 
introduced in the work of Mack and Salam [5], Ferrara, Gatto, Grillo and Parisi [6-11] and 
Polyakov [12]. In addition, we will need explicit expressions for the conformal blocks. Here 
we will be able to rely on the recent work of Dolan and Osborn [13-15]. 

While most of these ingredients were understood many years ago, until recently it was 
not known how to put them together in order to extract dynamical information about CFTs. 
This important know-how was developed in a series of recent papers [16-23]. That work 
was largely motivated by particle physics (in particular the theory of electroweak symmetry 
breaking) and concerned CFT in D = 4. However, the time is now ripe to transfer these 
techniques to D = 3. The cases D = 3 and D = 4 are similar in that the conformal algebra 
has finitely many generators (unlike in D = 2 where it has an infinite-dimensional extension, 
the Virasoro algebra). 

This paper is structured as follows. In Section 2 we review what is known about the 
operator content of the 3D Ising model. In Section 3 we discuss the conformal bootstrap 
approach to studying 3D CFTs, and in Section 4 we present an efficient method for com- 



puting the conformal partial waves appearing in four-point functions of scalars for CFTs in 
any dimension (including D = 3). In Section 5 we present bounds on 3D CFTs that follow 
from crossing symmetry and compare them to what is known about the 3D Ising model. 
Finally, we discuss our results and future directions for this program in Section 6. 



2 Operator Content of the 3D Ising Model 

We assume that the reader is familiar with the basic facts about the Ising model and the 
critical phenomena in general, see [24-27, 1]. 

In this paper, we will be aiming for a solution of the 3D Ising model in the continuum 
limit and at the critical temperature T = T^. While the 2D Ising model was solved exactly 
on the lattice and for any temperature by Onsager and Kaufman in the 1940 's, the 3D lattice 
case has resisted all attempts for an exact solution. Istrail [28] proved in 2000 that solving 
the 3D Ising model on the lattice is an NP-complete problem. However, this theorem does 
not exclude the possibility of finding a solution in the continuum limit. 

The standard way to think about the continuum theory is in terms of local operators (or 
fields). At T = T^, the theory has scale (and, as we discuss below, conformal) invariance, 
and each operator is characterized by its scaling dimension A and 0(3) spin. The operators 
of spin higher than 1 are traceless symmetric tensors. 

In Table 1 we list a few notable local operators, which split into odd and even sectors 
under the global Z2 symmetry (the Ising spin flip). The operators a and e are the lowest 
dimension Z2-odd and even scalars respectively — these are the continuum space versions of 
the Ising spin and of the product of two neighboring spins on the lattice. The two next- 
to-lowest scalars in each Z2-sector are called a' and e' . Their dimensions are related to 
the irrelevant critical exponents ua and u measuring corrections to scaling. The operator 
e" is analogously related to the next-to- leading Z2-even irrelevant exponent uj2- The stress 
tensor T^y has spin 2 and, as a consequence of being conserved, canonical dimension At = 3. 
The lowest-dimension spin 4 operator C^^kX has a small anomalous dimension, related to 
the critical exponent wnr measuring effects of rotational symmetry breaking on the cubic 
lattice. 



Operator 


Spin / 


Z2 


A 


Exponent 


a 





— 


0.5182(3) 


A = 1/2 + 77/2 


a' 





— 


>4.5 


A = 3 + wa 


e 





+ 


1.413(1) 


A = 3-l/z/ 


e' 





+ 


3.84(4) 


A = 3 + a; 


e" 





+ 


4.67(11) 


A = 3 + a;2 


T 


2 


+ 


3 


n/a 


^ IiukX 


4 


+ 


5.0208(12) 


A = 3 + wnr 



Table 1: Notable low- lying operators of the 3D Ising model at criticality. 



The approximate values of operator dimensions given in the table have been determined 
from a variety of theoretical techniques, most notably the e-expansion, high temperature 
expansion, and Monte-Carlo simulations; see p. 47 of Ref. [1] for a summary. The achieved 
precision is rather impressive for the lowest operator in each class, but quickly gets worse 
for the higher fields. While ultimately we would like to beat the old methods, it would be 
unwise to completely dismiss this known information and restart from scratch. Rather, we 
will be using it for guidance while sharpening our own methods. 

Among the old techniques, the e-expansion of Wilson and Fisher [2] deserves a separate 
comment. The well-known idea of this approach is that the 3D Ising critical point and the 
4D free scalar theory can be connected by a line of fixed points by allowing the dimension 
of space to vary continuously between 3 and 4. For D = 4 — e, the Wilson- Fisher fixed point 
is weakly coupled and the dimensions of local operators can be expanded order-by-order in 
e. For the most important operators, like a and e, these expansions have been extended to 
terms of order as high as e' [26] , requiring a five-loop perturbative field theory computation. 
However, as often happens in perturbation theory, the resulting series are only asymptotic. 
For the physically interesting case e = 1, their divergent nature already starts to show 
after the first couple of terms. Nevertheless, after appropriate resummation the e-expansion 
produces results in agreement with the other methods. So its basic hypothesis must be 
right, and can give useful qualitative information about the 3D Ising operator spectrum, 
even where accurate quantitative computations are missing. 

It is now time to bring up the conformal invariance of the critical point, conjectured 
by Polyakov [3]. This symmetry is left unused in the RG calculations leading to the e- 
expansion, and in most other existing techniques.^ This is because it only emerges at the 
critical point; it's not present along the flow. Conformal invariance seems to be a generic 
feature of criticality, but why exactly is not fully understood [31]. Recently there has been 
a renewed interest in the question of whether there exist interesting scale invariant but not 
conformal systems [32-37]. We will simply assume as a working hypothesis that the 3D 
Ising critical point is conformal. 

A nice experimental test of conformal invariance would be to measure the three-point 
function {cT{x)a{y)e{z)) on the lattice, to see if its functional form agrees with the one fixed 
by conformal symmetry [3] . We do not know if this has been done. 

Using 3D conformal invariance, local operators can be classified into primaries and 
descendants [5] . The primaries^ transform homogeneously under the finite-dimensional con- 
formal group, while the descendants are derivatives of primaries and transform accordingly. 
All operators listed in Table 1 are primaries. This is obvious for a and e — the lowest 
dimension scalars in each Z2-symmetry class. That a', e', e", C^^^x are all primaries and 
not derivative operators follows from the fact that they are associated with corrections to 
scaling, while adding a derivative operator to the Lagrangian has no effect. Finally, the 
stress tensor is always a primary. 

It can be seen that all operators in Table 1 have non-negative anomalous dimensions 

^Conformal invariance has been used in studies of critical 0{N) models in the large N limit [29, 30]. 
^These are usually called quasi-primaries in 2D CFTs. 



(by which we mean the difference between the operator dimension and the dimension of 
the lowest 3D free scalar theory operator with the same quantum numbers). This is not 
accidental, but is related to reflection positivity, which is the Euclidean space version of 
unitarity. Primaries in reflection positive (or unitary) CFTs are known to have non-negative 
anomalous dimensions [11, 38-41]: 

A>L)/2-l (/ = 0), A>l + D-2 (/>1). (2.1) 

The 3D Ising model is reflection positive on the lattice [42], and this property is inherited 
in the continuum limit, so that the 'unitarity bounds' (2.1) are respected. 

Can conformal symmetry be used to determine the local operator dimensions rather 
than to interpret the results obtained via other techniques? In 2D this was done long ago 
[4] using the Virasoro algebra. This also justified post factum the assumption of conformal 
invariance, since the critical exponents and other quantities agreed with the exact lattice 
solution. The Virasoro algebra does not extend to 3D, but in the next section we will 
describe a method which is applicable for any D. 



3 Conformal Bootstrap 

Primary operators in a CFT form an algebra under the Operator Product Expansion (OPE). 
This means that the product of two primary operators at nearby points can be replaced 
inside a correlation function by a series in other local operators times coordinate-dependent 
coefficient functions. Schematically, the OPE of two primaries has the form: 

0j(Xi)0j(x2) = ^ fijk C{Xi - X2, d2)(t)k{x2)- (3.1) 

k 

The differential operators C are fixed by conformal invariance, and only primary operators 
need to be included in the sum on the RHS. Here we are suppressing indices for clarity. In 
general, scalar operators as well as operators of nonzero spin will appear on the RHS. Fairly 
explicit expressions for the C's have been known since the 70's, at least in the case when 
the (pij are scalars and (f)k is a traceless symmetric tensor of arbitrary rank [6, 13], but we 
will not need them here. 

The numerical coefficients fijk are called structure constants, or OPE coefficients. These 
numbers, along with the dimensions and spins of all primary fields, comprise the 'CFT data' 
characterizing the algebra of local operators. 

The conformal bootstrap condition [8, 12, 4], shown schematically in Fig. 1, says that 
the operator algebra must be associative. In that figure we consider the correlator of four 
primaries 

{(t)l{xi)(t)2{x2)4>z{xz)<pA{XA)) (3.2) 

and use the OPE in the (12) (34) or (14)(23)-channel to reduce it to a sum of two-point 
functions. The answer should be the same, which gives a quadratic condition on the 



structure constants of the schematic form 

/ , /l2fc/34fe(- • •) = 2^ fl4kf23k{- ■ •) • (3.3) 

k k 

The (. . .) factors are functions of coordinates Xj, called conformal partial waves. They are 
produced by acting on the two-point function of the exchanged primary field 0^ with the 
differential operators C appearing in the OPE of two external primaries. Thus, they are also 
fixed by conformal invariance in terms of the dimensions and spins of the involved fields. 





4>:: 



Figure 1: The conformal bootstrap condition = associativity of the operator algebra. 



The dream of the conformal bootstrap is that the condition (3.3), when imposed on four- 
point functions of sufficiently many (all?) primary fields, should allow one to determine the 
CFT data and thus solve the CFT. Of course, there are presumably many different CFTs, 
and so one can expect some (discrete?) set of solutions. One of the criteria which will help 
us to select the solution representing the 3D Ising model is the global symmetry group, 
which must be Z2. 

Our method of dealing with the conformal bootstrap will require explicit knowledge of 
the conformal partial waves. In the next section we will gather the needed results. 



4 Conformal Blocks 

In this paper we will be imposing the bootstrap condition only on four-point functions of 
scalars. Conformal partial waves for such correlators were introduced in [7] and further 
studied in [9, 10]; they were also discussed in [12]. Recently, new deep results about them 
were obtained in [13-15]. Significant progress in understanding non-scalar conformal partial 
waves was made recently in [43] (building on [44]), which also contains a concise introduction 
to the concept. Below we'll normalize the scalar conformal partial waves as in [15]; see 
Appendix A for further details on our conventions. 

Consider a correlation function of four scalar primaries 0, of dimension A,, which is fixed 
by conformal invariance to have the form [3] 



Wiy l)V2y 2)V3y VVAK A)l l ^2 J \^^2j (^2^)i(A,+A,)^^2j|(A3+A4) ' ^ > 



where Xij = Xi — Xj, A^ = Aj — Aj, and g{u,v) is a function of the conformally invariant 
cross-ratios 

2 2 2 2 

^=^12^, t;=^M^. (4.2) 

X]^3X24 "''13"''24 

The conformal partial wave expansion in the (12) (34) channel gives a series representation 
for this function: 

g(u, V) = ^ fl2of-UO GA,liu, v) , (4.3) 

o 
where the sum is over the exchanged primaries O of dimension A and spin / and the functions 
Ga,i{u,v) are called conformal blocks. We must learn to compute them efficiently. 

In even dimensions, conformal blocks have relatively simple closed-form expressions in 
terms of hypergeometric functions [10, 13-15]. For example, the 2D and 4D blocks are given 

by: 

Gfj^u^v) = - [kA+i{z)kA-i{z) + {z^ z)] , 
G27\u,v) = T— ^^ [A:a+K^)^a-.-2(^) -{z^ -z)] , (4.4) 

i I X ^ /C' 

where 

k^{x) = x^/\F, (li^ - Au),lW + A,,); P;x) , (4.5) 

and the complex variable z and its complex conjugate z are related to u, v via 

u = zz, v = {l-z){l-z). (4.6) 

The meaning of the variable z is explained in Fig. 2. From the known analyticity properties 
of 2-^1, it follows that the conformal blocks are smooth single- valued functions in the z plane 
minus the origin and the (1, -|-oo) cut along the real axis. This is not accidental and should 
be valid for any D. By standard radial quantization reasoning (see [45], Sec. 2.9), the OPE 
by which the conformal blocks are defined is expected to converge as long as there is a sphere 
separating Xi and X2 from X3 and X4. This sphere degenerates into a plane and disappears 
precisely when z crosses the cut. 

We now pass to the results for general D, including the case D = 3 we are interested in, 
which are rather more complicated. From now on we consider only conformal blocks of four 
identical scalars, so that A12 = A34 = 0. In any dimension, such blocks depend only on 
the dimension and spin of the exchanged primary. For scalar exchange (/ = 0), conformal 
blocks have a double power series representation ([13], Eq. (2.32)): 

r (n, ,A - „^/2 \^ [(A/2),n,(A/2)n,+„] 

^0 m! n!(A + 1 - y)^(A)2™+n 

where (x)„ is the Pochhammer symbol. In this paper we will only use this representation 
aX z = z, in order to derive the closed form expression (4.10) given below. In principle the 
series converges absolutely in the region 

I -, < M < 1, ,, , 

- M, M > 1 , 




X4 ^ OO 




Figure 2: Using conformal freedom, three operators can be fixed at xi = 0, X3 = (1, 0, . . . , 0), 
X4 —7- 00, wliile the fourth point X2 can be assumed to he in the (12) plane. The variable z 
is then the complex coordinate of X2 in this plane, while z is its complex conjugate. Also 
shown: the conformal block analyticity cut (thick black), and the boundary of the absolute 
convergence region of the power series representation (4.7) (thin red). 



whose boundary is traced in red in Fig. 2. 

For exchanged operators of nonzero spin, the conformal blocks can be computed via vari- 
ous recursion relations. Some recursion relations previously appeared in [13], Eqs. (2.30), but 
these will not be useful for us since they express the blocks with equal external dimensions 
in terms of blocks where the external dimensions differ by an integer. 

In Appendix A, we exhibit a recursion relation which follows from the results of Ref. [15] 
and does not require shifts in the external dimensions. In general, this recursion involves 
taking derivatives of Ga,i{z,z), which is not very easy to perform numerically. However, 
along the line z = z the terms involving derivatives drop out and the recursion relation for 
the conformal block G/\^i{z) = G^^i{z,z) becomes extremely simple: 



{1 + D-3){2A + 2-D)GaAz) 



{D-2){A + 1-1)Ga. 



1-2 



z) + 



2- z 
~2^ 



(2/ + D-4)(A-D + 2)Ga+m-i(^) 



A(2/ + £) - 4)(A + 2 - D)(A + 3 - D)(A - / - D + 4)^ 



16(A + l-f)(A-f 



2)(/ - A + D - 5)(/ - A + £) - 3) 



G'A+2,i-2(^). (4.9) 



This recursion relation can easily compute all conformal blocks along the z = z line in 
terms of spin and 1 blocks. '^ On the other hand, as shown in Appendix B, the spin and 1 
blocks along the z = z line can be simply expressed in terms of generalized hypergeometric 



■^This works for general D. In D 
latter equality follows from (A. 3). 



3, one can instead recurse from Ga.o ^-nd Ga.-i = Ga.oi where the 



functions {a = D/2 — 1): 



2 \ A/2 X o 



Ga,o W = ^^ J 3F2 (^f , f , f - «; ^, A - «; ^^-^J , (4.10) 

A+l 

Ga,i(.) = ^ (^) ^ 3F2(^,^,^-«;f + l,A-a;^^). (4.11) 

These explicit expressions, together with the recursion relation (4.9), solve the problem 
of finding conformal blocks along the z = z line. What about z ^ z7 We should explain that 
in our numerical implementation of conformal bootstrap we will not actually use the values 
of conformal blocks at generic values of z. Instead, we will Taylor-expand the conformal 
bootstrap condition around the point z = z = 1/2. This is an approach which proved 
efficient in prior work in 4D and 2D, and we will pursue it here as well. So, we will have to 
evaluate derivatives of conformal blocks at the point z = z = 1/2, both along and transverse 
to the z = z line. 

Now, derivatives along the z = z line will be evaluated as follows. For the spin and 1 
conformal blocks we can take advantage of the fact that the 3F2 hypergeometric functions 
satisfy a third-order differential equation: 

(xiDa^Da2Da-i - I^oAi-i A2-1 j 8^^2(01, ^2, ^s; &1, h\ x) = , (4.12) 

where Dc = xdx + c. This equation can be used to obtain recursion relations which express 
the third-order and higher derivatives of the spin and 1 blocks in terms of their first and 
second derivatives. The values of the latter derivatives at z = z = 1/2 will be tabulated as 
a function of A. 

Derivatives of the higher spin blocks are then computed using the recursion relations 
following from (4.9). This completely settles the question of obtaining derivatives along the 
z = z line. 

In order to obtain the derivatives transverse to the z = z line, we'll take advantage of the 
fact that conformal partial waves are eigenfunctions of the quadratic Casimir operator of 
the conformal group, which implies that conformal blocks satisfy a second-order differential 
equation [14]: 

VGaAz, z) = \Ca^iGaj{z, z) , (4.13) 

where Ca,; = A(A -D) + l{l + D-2) and 

„, zzil — z) 



V = (l- z)z'd: 



2q2 



z'-iD 



d, + {z^z). (4.14) 



z — z 
Let us now make a change of variables: 

z = {a + Vb)/2, z={a-Vb)/2. (4.15) 

9 



The point z = z = 1/2 which interests us corresponds to a = 1, 6 = 0. Moreover, since 
conformal blocks are symmetric in ^ -H- ^, their power series expansion away from the z = z 
hne will contain only even powers of {^z — z\ and hence integer powers of h. In the new 
variables the differential operator T> takes the form 



V 



{2 - a]a-^ 



1 



{D-l)db + hd, 



+ {{2- 3a)b'dt + 



-{D - 9)a - a(3a - A)da -D + 3 



+ { -b^dadb + 



\{D-A)d. + \{2-3a)dl 



bdb-^Da^da+^{2-a)a^dl 



(4.16) 



where the terms have been grouped into lines according to how they change the power of b. 
The first line contains the leading terms, which lower the power of b by one unit. Notice that 
the leading terms generate a nonvanishing coefficient when acting on any positive power of 
b, as long as < a < 2 (which corresponds to < ^ < 1). Thus, in a neighborhood of this 
interval the Casimir differential equation (4.13) can be solved a la Cauchy-Kovalevskaya, 
recursively in a power series expansion in b using the known conformal blocks at 6 = as a 
boundary value. 

Let us denote the d^d]^ derivative of the conformal block evaluated at z = ^ = 1/2 by 
hm,n- Since we know the conformal blocks along the z = z line, we can compute all the 
derivatives hm,o- On the other hand, the Cauchy-Kovalevskaya argument above implies that 
there will be a recursion relation for hm,,n (with n > 0) in terms of /i^.n with lower values 
of n. The recursion relation is given in Appendix C and has the general structure: 



k 



m,n 



^ m{...)hm',n+ ^ [{...)hm',n-i + {n-l){...)h 



m'n—2\ 



(4.17) 



m'<m—l 



m'<m+2 



The appearance of m' up to m + 2 is related to the fact that, the Casimir equation being 
second-order, derivatives of up to second-order in a appear in the RHS of (4.16). The first 
term being proportional to m ensures that h^n terms generated by repeatedly applying the 
recursion are eventually reduced to hi^ and /io,o- This recursion then solves the problem of 
computing the conformal block derivatives transverse to the line z = z. 



5 Bounds and Consequences for the 3D Ising Model 

In this section we will use the bootstrap equations discussed above in order to derive rigorous 
bounds on 3D CFTs. When comparing these bounds to the 3D Ising model, we'll focus on 
constraints coming from the four-point function of the Ising spin operator (aacra). The 
conformal block expansion of this four-point function has the form 



g[u,v) 



^Pa,iG, 



AJlM,^^j 



PA,l = fli>0 



(5.1) 



where the sum runs over the dimensions and spins of all primary operators appearing in the 
a X a OPE. This OPE contains all of the Z2-even operators listed in Table 1, in addition 



10 



to infinitely many other even-spin operators. Note that odd-spin operators cannot appear 
because of Bose symmetry. The coefficients pa,/ appearing in the conformal block expansion 
are squares of the OPE coefficients, and are thus constrained to be positive. 

The conformal bootstrap equation (3.3) takes a particularly simple form for this corre- 
lator, since the (12) (34) and (14) (23) channel involve the same OPE coefficients. It can be 
stated as a crossing symmetry constraint on the function g(u,v): 

v^^g^u, v) = u^^giv, u) . (5.2) 

Substituting the conformal block decomposition, we get an equation 

^A. _^A. ^ Y^p^, [v'^^G^^iiu.v) -u'^^G^^iiv^u)] , (5.3) 

where ^ is the sum over all operators except the unit operator, whose contribution has 
been separated in the LHS. It was shown in [16] and confirmed in subsequent work [17-23] 
that this type of equation can be used to extract dynamical information about 4D and 2D 
CFTs. We will now apply the same methods in 3D. 

First, we will Taylor-expand (5.3) around the point z = z = 1/2 up to some large 
fixed order. That this is a reasonable point to expand around follows from the fact that 
it is democratic with respect to the direct and crossed channels in the conformal block 
decomposition: by making a conformal transformation the four points can be put at the 
vertices of a square. 

The Taylor-expanded (5.3) can be viewed as a finite system of linear equations (one for 
each Taylor coefficient) for a large (strictly speaking infinite) number of variables pA,i- A 
priori, there is one variable pa,^ for each pair (A, /) consistent with the unitarity bounds 
(2.1). However, one may wish to posit additional constraints on the spectrum (such as 
assumptions about gaps). Below we will study which of these constraints are consistent 
with the existence of a solution. 

This system of linear equations should also be augmented by inequalities expressing the 
fact that variables pa,« are non-negative. Fortunately, problems involving linear inequalities 
are almost as tractable as pure systems of linear equalities. These problems form a chapter 
of linear algebra called linear programming, and there exist efficient algorithms for solving 
them (such as Dantzig's simplex method or interior point methods). Once the additional 
constraints on the spectrum are specified, one can use linear programming methods to find 
out if the system has a solution. If the answer is negative, a CFT with such a spectrum 
cannot exist. The details of our implementation of this problem are given in Appendix D. 



5.1 Bounds on A 



£ 



We are now ready to start asking concrete questions about the 3D Ising CFT to which we 
can give unambiguous answers. The first question is as follows. Let's be agnostic about the 
dimension of the spin field, allowing it to vary in the interval 0.5 < A^. < 0.8. The lower 
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0.50 0.55 0.60 0.65 0.70 0.75 0.80 



Figure 3: Shaded: the part of the (Aq-, Ag) plane allowed by the crossing symmetry constraint 
(5.3). The boundary of this region has a kink remarkably close to the known 3D Ising model 
operator dimensions (the tip of the arrow) . The zoom of the dashed rectangle area is shown in 
Fig. 4. This plot was obtained with the algorithm described in Appendix D with rimax = 11- 



end of this interval is fixed by the unitarity bound, while the upper end has been chosen 
arbitrarily. For each Ao- in this range, we ask: What is the maximal A^ allowed by (5.3)? 

The result is plotted in Fig. 3: only the points (A^, A^) in the shaded region are allowed.^ 
Just like similar plots in 4D and 2D [16, 17, 23] the curve bounding the allowed region starts 
at the free theory point and rises steadily. Moreover, just like in 2D [17] the curve shows a 
kink whose position looks remarkably close to the Ising model point. ^ This is better seen in 
Fig. 4 where we zoom in on the kink region. The boundary of the allowed region intersects 
the red rectangle drawn using the A^ and A^ error bands given in Table 1. 



1.44 
1.43 
1.42 
1.41 
1.40 
1.39 



1.38^— 
0.510 




0.515 



0.520 



0.525 



0.530 



Figure 4: The zoom of the dashed rectangle area from Fig. 3. The small red rectangle is 
drawn using the Ao- and A^ error bands given in Table 1. 



From this comparison, we can draw two solid conclusions. First of all, the old results 
for the allowed dimensions are not inconsistent with conformal invariance, though they are 



■*To avoid possible confusion: we show only the upper boundary of the allowed region. 0.5 < A^ < 1 is 
also a priori allowed. 

^In contrast, the 4D dimension bounds do not show kinks, except in supersynimetric theories [23]. 
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based on completely different techniques. Second, we can rigorously rule out about half of 
the (Ao-, Ae) rectangle allowed by the table. It seems that the 3D Ising model lies remarkably 
close to the boundary of the allowed region, if not on the boundary. At present we don't 
have an explanation of why this had to be the case. 



5.2 Bounds Assuming a Gap Between £ and e' 

We will next give a series of plots showing the impact of assuming a gap in the Z2-even 
scalar spectrum (as proposed in [46]). In other words, we will impose that the first operator 
after e has dimension A^/ above a certain value. 

Going from weaker to stronger, we will consider three constraints: A^/ > 3,3.4,3.8. 
Thus, we will ask: What is the region of the (Ao-, A^) plane allowed by (5.3) when this extra 
constraint is taken into account? 

The weakest of the three assumptions, A^/ > 3, has been chosen since it can be justified 
experimentally: we know that the 3D Ising critical point is reached by fine-tuning just one 
parameter (the temperature). Therefore, it has just one relevant Z2-even scalar, £, while 
e',e" etc. must be irrelevant. As we see in Fig. 5(a), this piece of information allows to 
exclude a fair part of the region allowed by Fig. 3. Unfortunately, close to the 3D Ising we 
do not gain constraining power: the new and the old bounds coincide there. 

On the other hand, the stronger assumptions A^/ > 3.4, 3.8 exclude a much larger portion 
of dimension space, carving out an allowed region with two branches; see Figs. 5(b,c). The 
upper branch seems to end at the 3D Ising point, while the lower branch terminates near the 
free theory. It is simple to understand why the intermediate region should not be allowed 
- assuming a gap A^/ > A^, should exclude the gaussian line A^ = 2Ao- up to a dimension 
of Ao- = A*/2 — 1, since the spectrum of this solution is 2Ao- + 2n + / for integer n. Our 
bounds are slightly weaker than that. 

Zooming in on the tip near the 3D Ising point, we see that the allowed region in Fig. 5(c) 
barely intersects with the red rectangle. Were we to assume even larger gaps, the intersection 
would eventually disappear altogether. We performed this analysis and found that this 
happens for A^/ > 3.840(2). This result rules out the upper half of the A^/ range allowed by 
Table 1, assuming that the more accurate determinations of A^. and A^ in the same table 
are correct. 

The same phenomenon is seen in a slightly different way in Fig. 6. Here we compute the 
maximal allowed A^/ under the condition that A^ has already been fixed to the maximal 
value allowed by Fig. 3. Notice the rapid growth of the A^/ bound just below the 3D Ising 
model a dimension, which allows e' to become irrelevant. Similar growth has been observed 
in the 2D case in [46]. Around the 3D Ising A^- the bound is A^/ < 3.84, consistent with 
the value cited above. 

This story illustrates how the conformal bootstrap equation imposes nontrivial depen- 
dencies between various operator dimensions. Once some dimensions are determined, the 
other ones are no longer arbitrary. Such interrelations are probably not easy to see from 
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Allowed Region Assuming A{e')>3 



(Zoomed) Allowed Region Assuming A(e')>3 
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Allowed Region Assuming A(e')£3.4 



(Zoomed) Allowed Region Assuming A(e')>3.4 
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Allowed Region Assuming A(e')£3.8 
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(b) 



(Zoomed) Allowed Region Assuming A(e')>3.8 
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Figure 5: Same as Figs. 3, 4, but imposing the extra constraints A^/ > {3,3.4,3.8}. 
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Figure 6: The bound on A^/ under the condition that A^ has already been fixed to the 
maximal value allowed by Fig. 3. Here nmax = 10 (see Appendix D). The width of the vertical 
red line marking the 3D Ising value of Aq- is about five times the error band in Table 1. 



the renormalization group point of view. For instance, when using the e-expansion, each of 
the operator dimensions hsted in Table 1 requires an independent computation. 



5.3 Bounds on the Gap in the Spin 2 Sector 



The above discussion concerned the scalar sector of the 3D Ising model, but eventually we 
would like to also constrain operators with nonvanishing spin. For a first try, let's study 
here the gap in the spin 2 sector. The first spin 2 operator in the a x cr OPE is the stress 
tensor T^^^, and we will be interested in the dimension of the second one, call it T' . 

In Fig. 7 we give a rigorous upper bound on At> following from the crossing symmetry 
constraint (5.3). The bound is shown as a function of Ao- only, and is in this sense analogous 
to our first bound in Fig. 3. Unlike for the case of e' studied in the previous section, we 
found that the bound on T' is only very weakly correlated with the value of Ag, and so we 
do not show separately the allowed regions in the (Ao-, A^) plane. 

The At' bound shows fascinating behavior which is the opposite to that of Fig. 6. It 
has a plateau at A^' ~ 5.7 for low A^- and suddenly drops to much lower values Aj" ~ 3.5 
as the dimension of a is increased. To begin with, this implies that any moderate gap in 
the T' dimension, e.g. A^-/ > 4, leads to a sharp upper bound on A„. Taken together with 
the plots in Fig. 5, one then obtains very small closed regions in the (A^, A^) plane. 

Moreover, the sudden drop in the At> bound happens precisely when Ao- passes the 3D 
Ising value. The actual bound there is: 



A. 



0.518 



At' < 5.6 



(5.4) 



Unfortunately, Table 1 is mute about Ay/ as we are not aware of any prior studies. However, 
we can get a rough estimate of this dimension by interpolating between 2D and 4D. In the 
4D free scalar theory the first Z2-even spin 2 operator after the stress tensor is 



r 



'T, 



flU • 



(4D) 



(5.5) 
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Figure 7: Upper bound on the dimension of the second spin 2 operator T' from the crossing 
symmetry constraint (5.3). The algoritlim from Appendix D was used with n-max = 10. The 
3D Ising vertical red hne is five times wider than the error band in Table 1. We do not show 
the region of A^ close to the unitarity bound, which is subject to numerical instabilities. 



which has dimension 6. To be more precise, in the free scalar theory this operator is 
decoupled from the 0x0 OPE, but we expect it to couple in the Wilson- Fischer fixed point 
in 4 — e dimensions. 



In the 2D Ising model the first such quasiprimary operator is 



r = (L_ 



3-^-2)-^-2l 



(2D) 



(5.6) 



lL'_,)e 



again of dimension 6. Notice that another 2D candidate spin 2 quasiprimary, (-L_2 

of dimension 3, is a null state since the field e = 02, i is degenerate on level 2 in the 2D Ising 

model. 

Assuming as usual that the 2D Ising and the 4D free scalar are continuously connected 
by the line of Wilson-Fischer fixed points to which the 3D Ising model also belongs, we 
expect by interpolation that Aj" ~ 6 in 3D, not far from the upper end of the range allowed 
by the rigorous bound (5.4). 



5.4 Bounds on Higher Spin Primaries 

In addition to bounding operators in the scalar and spin 2 sectors, we can also attempt to 
place bounds on higher spin primaries in the a x a OPE. The first such operator in the 3D 
Ising model is the spin 4 operator C^^kX- This operator is interesting because it controls 
the leading effects of rotational symmetry breaking when the 3D Ising model is placed on a 
cubic lattice. The corresponding perturbation of the CFT Lagrangian can be written as 



^-CcFT OC Ciiii + C2222 + Q 



3333 



(5.7) 



Because of this connection with phenomenology, the dimension of C has been computed 
rather precisely: Ac ^ 5.0208(12) ([47], Eq. (4.9)). 
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Figure 8: Upper bound on the dimension of tiie first spin 4 operator in the a x a OPE 
from the crossing symmetry constraint (5.3). The algorithm from Appendix D was used with 
'^max = 10. The tip of the arrow shows the point (Aq-, Ac) with the 3D Ising model values 
from Table 1. The dashed line is the gaussian solution A4 = 2Ao- + 4. 



In Fig. 8 we give a rigorous upper bound on Ac following from crossing symmetry and 
unitarity, making no other assumptions about the spectrum. While this bound passes above 
the value of Ac in the 3D Ising model, this is easily understood by the fact that the gaussian 
solution to crossing symmetry has Ac = 2Ao- + 4, which must be respected by our bound. 
The interesting and highly nontrivial statement is then that the gaussian solution seems to 
essentially saturate the bound. The bound that we find is fit well by the curve: 



A max 



(2A, + 4) + 0.1176 (A, - 1/2)' + O ((A, - 1/2)^) 



(5.8) 



so that one can see that linear deviations are not required, and quadratic deviations are at 
least somewhat suppressed. It is tempting to conjecture that the optimal bound (taking 
'^max -> 00) will exactly follow the gaussian line. It will be important in future studies 
to closely examine behavior of the bound at even larger external dimensions, to better 
understand whether deviations from this conjectured behavior are allowed. 

Does this behavior of closely following the gaussian line hold for higher spins? To explore 
this, in Fig. 9 we show the analogous upper bound on the lowest- dimension spin 6 operator 
in the a x a OPE. This operator would control breaking of rotational symmetry on the 
tetrahedral lattice, but we are not aware of prior 3D studies of its dimension. Again we see 
that the bound closely follows the gaussian line Ag = 2Ao- + 6, with a fit: 



(2A, + 6) + 0.1307 (A, - 1/2)' + O ((A, - 1/2)') 



(5.9) 



so that again quadratic deviations are suppressed. We have verified that this trend continues 
for operators of spin 8 and 10. 

An important feature of these bounds is that they approach the dimensions of spin / 
conserved currents A; = / + 1 as Aq- — > 1/2. It is well known that theories of free scalars 
contain higher spin conserved currents. Our bound shows that theories containing almost- 
free scalars necessarily contain higher spin operators that are almost conserved currents. 
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Figure 9: Upper bound on the dimension of the first spin 6 operator in the a x a OPE 
from the crossing symmetry constraint (5.3). The algorithm from Appendix D was used with 
JT-max = 10. The dashed hne is the gaussian solution Ag = 2Ao- + 6. 



A CFT version of the Coleman-Mandula theorem proved recently in [48] shows that theories 
containing higher spin currents and a finite central charge necessarily have the correlation 
functions of free field operators. This implies that we should also be able to derive a lower 
bound on the dimensions of higher spin operators, perhaps under the assumption of a finite 
central charge. It would be also interesting to connect these studies with an old result of 
Nachtmann [49] that in a unitary theory the leading twists 



Ti = Ai-{l + D-2), 



(5.10) 



where A; is the dimension of the lowest spin / operator, must form a nondecreasing and 
convex upward sequence for / > 2. We leave exploration of these very interesting directions 
to future work. 

These bounds are also particularly interesting in the context of the AdS/CFT correspon- 
dence, since they place tight constraints on 0(1/A^^) corrections to the dimensions of double- 
trace operators. Concretely, free scalars in AdS give rise to spin-/ double-trace operators 
with gaussian dimensions 2Ao-+2n+Z for integer n, while bulk interactions generate 0(1/A^^) 
corrections to these dimensions. Some explicit examples of these corrections were studied, 
e.g., in [50, 51]. If our conjecture that the gaussian solution saturates the bound is true, then 
the bounds forbid bulk interactions that generate positive corrections to these dimensions, 
which in turn may imply positivity constraints on (higher derivative) interactions in AdS. 
Such constraints could then be related to the constraints on higher derivative interactions 
studied in [52]. This is clearly another direction worth studying in future work. 

Finally, let us mention that similar bounds can be derived on the lowest dimension spin 
2 operator in (non-local) theories where a stress tensor does not appear in the a x a OPE. 
This bound (Fig. 10) shows similar features to the higher spin bounds. 

Such non-local theories may be interesting for several reasons. First, they commonly 
arise in statistical mechanics as models of long-range critical behavior. One much studied 
example is the critical point of the long-range Ising model, defined by a lattice Hamiltonian 
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Figure 10: Upper bound on the dimension of the first spin 2 operator in the a x a OPE 
from the crossing symmetry constraint (5.3) in non-local theories without a stress tensor. The 
algorithm from Appendix D was used with rimax = 10. The dashed line is the gaussian solution 
A2 = 2A^ + 2. 



with a power-law spin-spin interaction: 

n 



E 



SiSj 



(5.11) 



i,j V 



The precise universality class of this model depends on the value of 7. According to classic 
results [53, 54] supported by Monte-Carlo simulations [55]^ , there are three regions. For 
7 sufficiently small, namely 7 < D/2, the critical point is the gaussian model with the 
spin-field dimension determined by the naive continuous limit of (5.11): A^- = {D — 7)/2. 
Then there is an intermediate region, and finally the region of large 7, in which the model 
belongs to the usual, short-range, Ising model universality class and the critical exponents 
do not depend on 7. The boundary between the intermediate and short-range region lies at 
7 = D — 2Ajf™^, determined by the short-range Ising model spin-field dimension. This can 
be also understood by studying stability of the short-range Ising model with respect to non- 
local perturbations. Analogously, the boundary between the gaussian and the intermediate 
region lies at the value of 7 for which the operator o"^ becomes marginal. 

In the intermediate region, the a dimension is still given by the mean-field formula A^. = 
(D — 7)/2, but the dimensions of other operators, such as e, have nontrivial dependences on 
7 deviating from the gaussian values. So these fixed points are interacting. Because of their 
origin as relevant perturbations of the non-local gaussian scalar theory, they are expected 
to have conformal symmetry (and not just scale invariance), but not a stress-tensor. It is 
for such non-local CFTs that our bound in Fig. 10 may be of interest. 

Another reason to be interested in theories without a stress tensor is that they realize 
a simpler case of AdS/CFT, in which bulk gravity is decoupled, so that the AdS metric is 

^Note added: very recent accurate Monte-Carlo simulations [56] indicate deviations near the boundary 
of the intermediate and the short-range regions from the estabhshed picture as described below. 



19 



viewed as a fixed non-fluctuating background/ This may be useful wlien one is interested 
in aspects of tlie correspondence which are not necessarily related to gravity, as e.g. in 
[50]. Also, removing gravity allows one to flnd nontrivial UV-complete AdS/CFT examples 
which are purely fleld-theoretic (no strings): any UV-complete quantum fleld theory on 
the AdS/j+i background can be interpreted as providing a dual description to a non-local 
£)- dimensional CFT on the boundary. 



5.5 Bounds on the Central Charge 

Our final application concerns the central charge Ct of the 3D Ising model, defined for an 
arbitrary D as the coefficient of the canonically normalized stress tensor two-point function: 



Ct 1 



i^yJ^fiX-^ua I -i fiaJ^uX) r^'^ fiv'^Xa 



V = V - 2x^x,/a;2 , So = 2n^/yT{D/2) . (5. 12) 

It seems that the 3D Ising central charge has been computed only to the second order in 
the e-expansion, with the result [58-60, 30] 

CT/C'r = ^-^e' + 0{e'), (5.13) 

where C^'^'^ = D/{D — 1) is the free scalar field central charge. Substituting e — )■ 1 and 
neglecting the unknown higher-order terms, this estimate would suggest that Ct/C^^^ is 
very close to 1, around 0.98 or so. 

In our method, we can get control over Ct because the stress tensor conformal block 
enters the crossing symmetry constraint (5.3) with a Ct dependent coefficient:^ 

D A2 

Following [19, 20], the conformal bootstrap can be used to bound the coefficient P3.2 from 
above, which bounds the central charge from below. In Fig. 11 we show the lower bound on 
Ct as a function of A^.. We see that the bound has a distinctive minimum close to the 3D 
Ising value of the a dimension. The position of the minimum corresponds to Ct/C^'^^ ~ 0.94. 
One may also redo the plot in Fig. 11 making some assumption about A^, like that A^ > Ao-. 
The most aggressive assumption would be to fix A^ to the maximal value allowed by the 
upper bound in Fig. 3. One finds that the shape of the bound on Ct is very weakly 
dependent on these assumptions, but that the minimum moves to the right, even closer to 
the 3D Ising A^,, and slightly higher up to Ct/C^'^*' ^ 0.95. 

We believe that the observed minimum in the Ct lower bound is not accidental, but 
must be close to the true value of Ct-^ This would imply a small but noticeable discrepancy 



^Such theories may alternately be viewed as the starting point for "constructive holography" by defining 
a CFT perturbatively around the c — >■ 00 point as done in e.g. [57]. 

^The prefactor is different from [20] due to the different conformal block normalization, see Eq. (A. 2). 

^In 2D, a similar analysis reproduces the exact value of the 2D Ising model central charge with 10"'' 
accuracy [61]. 
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with the e-expansion estimate of Ct-, which can be attributed to the unknown higher-order 
terms. In fact, we can also derive upper bounds on Ct in presence of a gap between T 
and T' . The strength of these bounds depends on the assumption about the gap, and for 
T' close to the maximal value allowed by (5.4) would rigorously rule out the e-expansion 
estimate. We leave full exploration of such upper bound bounds to future work. 




0.60 



Figure 11: The lower bound on Ct as a function of Ao-. The plot was obtained with 
?^max = 11. The 3D Ising vertical red line is five times wider than the error band in Table 1. 



6 Discussion 

The results of the previous section have many implications whose importance is hard to 
overestimate. First, all of our bounds are consistent with everything that was previously 
known about the critical exponents of the 3D Ising model, as computed via RG methods and 
measured in experiments and Monte Carlo simulations. We should take this a very strong 
evidence that the 3D Ising model has a full conformal symmetry, justifying post factum the 
use of conformal symmetry in studying this theory. It would be good to further test the 
conformal invariance experimentally or on the lattice, for example by measuring the form 
of the 3-point functions. One can also compare any new measurements (e.g., of the central 
charge) against the constraints obtained using the methods in this paper. 

It is worth emphasizing that the bootstrap approach to studying 3D CFTs taken in this 
paper has, in principle, a significant advantage over other methods - at every step in the 
program we can present constraints that are completely rigorous (up to numerical errors 
that can be made arbitrarily small). This is significantly better than the usual situation 
in field theory computations, where one computes the first several terms in a series (say, 
the e-expansion or a loop expansion) and one can only estimate the errors from neglecting 
higher terms. It is also an advantage over lattice simulations, where it can be very difficult 
to gain control over errors induced by discretizing the theory. 

In this paper we have imposed only the first and the simplest of the infinitely many 
bootstrap conditions - the one following from the crossing symmetry of the a four-point 
function. It turns out that this condition alone carves out a significant portion of the 
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operator dimension space. The 3D Ising model seems to lie on the boundary of the allowed 
region, and at a rather special point - a corner. This empirical fact suggests that the model is 
algebraically special, for two reasons. First, the crossing symmetry constraint is expected to 
allow fewer solutions at the boundary of the allowed region as compared to the bulk, perhaps 
just a unique solution. Second, the non-analytic behavior of the bound at a corner point 
can be attributed to rapid rearrangements of the operator spectrum [46]. Indeed, Figs. 6 
and 7 show rapid changes happening for the next-to-leading operator dimensions in the 
scalar and spin 2 sectors. Such spectrum rearrangements signal linear (near-) degeneracies 
among various conformal blocks. It is very important to explore this phenomenon in detail 
as it offers tantalizing hope for distilling some analytical understanding of the 3D Ising 
model dimensions from our numerical approach. More generally, the fact that some special 
theories seem to lie at the edge of the region allowed by crossing symmetry may suggest a 
new classification scheme for understanding CFTs in D > 2. 

Furthermore, it is intriguing that most of our bounds (not just A^) seem to be essentially 
saturated by the values realized in the 3D Ising model. This fact suggests the strategy of 
determining the spectrum recursively: first fix e at the maximal allowed dimension, then e' 
at the maximal allowed dimension given A^, etc. We hope to explore the viability of this 
approach (perhaps also including gaps in higher spin operators) in future work. 

Another badly needed development is to add conformal bootstrap constraints coming 
from other correlators, which can lead to interesting interplay. For example, we would like 
to include (aeae) expanding in the axe channel, since this expansion will be crossing 
symmetric. Moreover, the conformal block of a will appear with the same coefficient /^^^ 
as the conformal block of e in the analysis of (aaaa). It is also interesting to include (eeee) 
whose expansion involves the same Z2-even operators as {aaaa). Due to (5.14), the stress 
tensor will appear in both expansions with related coefficients. One can also consider 4- 
point functions containing the stress tensor, where the recent results of [43] on conformal 
blocks for external operators with spin can be used. 

Another future task is to study 3D CFTs with larger global symmetry groups, such as 
0{N) symmetry. The general theory of analyzing bootstrap constraints in the presence of a 
continuous global symmetry was given in [21]. The equations look more difficult as the OPE 
contributions should be classified into various representations and the crossing symmetry 
transformation involves a Fierz matrix. Nevertheless, there are always as many equations 
as representation channels and the total constraining power is expected to be comparable 
to the Z2-symmetric case. In 4D, this has been convincingly demonstrated in [22, 23], 
where many strong bounds for 0{N) and SU{N) symmetric CFTs have been obtained. It 
would be interesting to generalize these methods to 3D and see how the resulting bounds 
on operator dimensions compare to what is known about the 0(A^)-vector models. 

Cross-fertilizing in the opposite direction, it is worth applying in 4D what we have learned 
in this paper in the 3D context - how interesting it is to study the effects of gaps in the 
operator spectrum. In addition, we should stress that the recursion relations for conformal 
blocks exhibited in this paper are valid for any space-time dimension D. Thus, we can use 
them to numerically compute conformal blocks in 4 — e dimensions for different values of 
e, where we can make contact with operator dimensions and OPE coefficients computed 
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perturbatively in the e-expansion. 

Our results and discussion in Section 5.4 show that one can also learn interesting state- 
ments about higher spin operators from crossing symmetry. It will be interesting to explore 
what can be learned further, particularly in the context of the AdS/CFT correspondence, 
where for example the 0(A^)-vector models in the large A^ limit are described by higher spin 
gauge theories in AdS4 [62, 63]. There is clearly still much to be learned about the role that 
higher spin operators play in ensuring consistency of the theory, and about how gaps in the 
lower-spin spectrum affect what these operators are allowed to do. 

Overall, our results fly in the face of the prevailing opinion that above two dimensions 
conformal symmetry by itself is not sufficiently restrictive to solve models. Clearly, the 
conformal bootstrap in D > 2 works. We have not yet solved the 3D Ising model, but we 
have definitely cornered it. 
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A Recursion Relations at Fixed External Dimensions 

Our conformal blocks are the same as the functions -Fai,A2 of R-^f. [15]: 

Ga,i = i^A,A„ Ai = i(A + /), A2 = ^(A - . (A.l) 

This normahzation is different from the one used in a number of previous works. E.g., in 
[13, 43] conformal blocks contain an extra factor of^° 

^ ^' a = 1. (A.2) 



(-2)H«)/ ' " 2 

This follows by comparing Eq. (2.25-29) of [15] with Eq. (2.22) of [43]. 

We note in passing one reason for using the new normalization: once conformal blocks 
are analytically continued to all real /, one has a symmetry relation ([15], Eq. (4.10)) 

Ga,i = Ga-i~d+2 ■ (A. 3) 

In particular, we have Ga.-i = G^a in 3D, which can be useful as explained in footnote 3. 

Below we consider only the case A12 = A34 = 0, which corresponds to setting a = b = 
in the notation of [15]. Denote 



/3p = 77;^ f' , ,w D,^z\l-z)d'i-z%, (A.4) 



J-o = - + --l, T, = il-z)d, + il-z)d,, jr^ = L^^D,-D,) . (A.5) 

z z zz 

It was shown in [15] that TiF\^\^ can be expressed as linear combinations of -Fa' a; • More 
specifically, we have (see [15], Eqs. (4.28), (4.29), (4.32)) 



/ + 2a 



-^0 F\i\2 — ~i~, FaiA2-i + 1—. — Fxi-iX2 



I (A-l)(A-2a) (l + 2a ^ ^ R p \ ( A 6) 

+ (A-l-«)(A-a) ^TT^ /^Aii^A.+iA, + Y^ /3a.-„Fa.A2+i J , (A.6) 

l + 2a I (A-l)(A-2a) 

•ri -TAiAa - -r, A2 -FA1A2-1 + -r, — (Ai + a) Fai_ia2 + -ri 1 rrr r x 

l + a l + a [A — 1 — a){A — a) 

X ( , , (-Ai + a + 1) /3aiFai+ia2 + ^ — (-A2 + 2a + 1) /3A2-ai^AiA2+i 
\ / + a l + a 

■(A.7) 
l{l + 2a) 



J-2Fa,a2 = (A-1 



/ + a 

(A-2a)(A- l-2a) 
(A-l-a)(A-a) 



FA1A2-I — -^Ai-1A2 



(/3ai -Fa-^+iA2 — (3x2-a FaiA2+i) 



(A.^ 



^^a was called e in [15], but we already have two other epsilons in this paper. 

24 



where A = Ai + A2, / = Ai — A2. 

Let us now view (A. 6) and (A. 8) as a linear 2x2 system for the spin / + 1 conformal 
blocks Fai+ia2 and FaiA2-i- Eliminating one of these, say Fa^+iaj, we get a recursion relation 
expressing the remaining spin / + 1 block in terms of spin / and spin / — 1 blocks only. Shifting 
the spin by one and passing to the Ga,/ notation, this relation takes the form: 

(A-«)(/ + 2«-l) a{^ + l-l) If 7,\ 



l + a-1 ' l + a-1 ' 2V /-I, 

A(A-2^)(A-2a + l) ' 

" (A-«)(A-« + l) /^|(A-^+2-2.)^A+2,-2 . (A.9) 

When we specialize to the line z = z, the term involving J-2 vanishes. We are then left with 
a nonderivative recursion relation, Eq. (4.9) of the main text. 

Alternatively, we can apply the same logic to the system formed by (A. 6) and (A. 7). 
Eliminating again Fx^+ix^ in favor of -FaiA2-1) shifting the spin by one and passing to the 
Ga,i notation, we get: 

(A - a)(/ + 2a - 1)^^^_ ^(l^^^l_2a- 2)7, + 7^ G.,„_. 



(/ + a-l) • V2 

- (^ - 1) ((A - «)(A -t + i)/5|(a-/+2-2.)G'a+2,-2 + 7^^3T^A,.-2j • (A.IO) 

Recursions (A.9) and (A.IO) have complementary advantages. The first one becomes 
nonderivative at 2; = ^ and can be used to compute high spin blocks on this line efficiently. 
However, it needs both / = and / = 1 blocks to start up (except in D = 3 where it can be 
started from / = and / = — 1, but we would like a framework which works in any D). On 
the other hand, recursion (A.IO) has spin / — 2 blocks entering with a factor (/ — 1) and can 
be started up with just spin 0. In Appendix B, we'll use (A.IO) to compute spin 1 blocks 
ai z = z from spin 0, but switch to the nonderivative recursion (A.9) for higher spins. 

B Scalar and Spin 1 Blocks at z = z 

In this appendix we'll derive formulas for the spin and 1 conformal blocks aX z = z for 
equal external dimensions. 

We start with the double series expansion (4.7) for the scalar conformal block. Perform- 
ing the summation in n, we get (a = D /2 — 1) 

°o nA\ w / A A \ 

^A,o = E UA^ (i ^«^+'"2FiU + -,m + -;2m + A;l-t; . (B.l) 

We now replace 2-^"! by its Euler integral representation 

T(c) r^ t^-^d-ty-^-^ 

^^'("' '• '■' ^> = m^) I * (1-4° • '^-'^ 
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The series in m under the integral sign turns out to be hypergeometric in the variable 

_ (1 - t)tu 

so that we find: 

Now let us use the hypergeometric identity 

2Fi(a, b; c; x) = (1 - x)"^ 2-^1 f c - a, 6; c; j . (B.5) 

The resulting expression factorizes nicely in terms of z and z: 

G.o = ^[' -^Y-f^ ,F, f ^, ^ - a; A - a; -y] , (B.6) 

^•° r(f)'yo t{i-t) ' 'V2' 2 ' ; ' ^ ^ 

y^^= '^'-'>'' . (B.7) 

Now replace 2-^1 by its defining power series expansion in {—Y) and integrate the series term 
by term. For z = z, the resulting integrals are of the form (4.7) and give hypergeometric 
functions 2-^1 (A + 2n, A/2 + n; A + 2n; z), which are elementary. ^-^ So we get: 

Ga,oU«- - (^^^ j 1. „!(A).„(A-a)„ [v^l ) ' <^-») 

Expressing (A)2n via the duplication formula for the F function, the series is recognized to 
be of the 3F2 type, and we get precisely Eq. (4.10). 

Is there a similar closed form representation for generic unequal external dimensions or, 
more specifically, for generic A12 = A34 7^ (as would be needed for the crossing symmetry 
analysis of the (aecre) correlator)? The following reasoning shows that this may be difficult. 
For D = 2, Eq. (4.10) can be derived starting from the explicit expression (4.4), passing to 
the variable z^/(4(z — 1)) via the identity 

2F,ia,b,2b,z) = (1 - I) (1 - ^)-^ 2F1 (^t?^, ^-f,b+l; ^) , (B.9) 

and then aiming for Clausen's formula ([64], Sec. 4.3) to express the square of a 2-^1 as a 
3F2. However, Eq. (B.9) is not useful for generic unequal dimensions. 

Passing to the spin 1 case, the idea is to use the second recursion relation (A. 10) which 
expresses spin 1 blocks via the spin ones. This relation can be restricted to the z = z line, 
as the differential operator J^i acts within the line; for Z = 1 it gives 



GaA^) 



2{A-a) 



^(A-2«-l) + (1-^)9, 

2z 



G'a+i,o(^). (B.IO) 



^^For z ^ z we would have obtained a series in Appel Fi functions. 
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Substituting the spin closed form expression (4.10), we find 



9 N ^±i 
vi \ 2 



gA.i(^) = ^^2-l)4 r^J [2/^. + A-«]/(2/), (B.ll) 

where f{y) (with y = z^/4:{z — 1)) is the 3F2 function entering the expression for Ga+i,o{z). 
Eq. (4.11) then follows, since the 3F2 function satisfies 

[ydy + 62 - i]3^2(ai, ^2, %; &i, ^2; y) = {h - 1)3^^2(^1, ^2, %; &i, &2 - 1; y) • (B.12) 



C Recursion Relation for the Transverse Derivatives 

The following recursion relation for hm,n can be derived by applying i9™9^ to the Casimir 
equation (4.13) written in the a, h coordinates, setting a — > 1, 6 — )■ 0, and shifting n — t- n — 1: 

2{D + 2n - 3)/i™,„ = 

2m{D + 2n- 3)[-hm-i,n + (m - l)hm-2,n + (m - l)(m - 2)hm-3,n] 

- hm+2,n-l + (D - m - An + 4:)hm+l,n-l 

+ [2Ca,/ + 2D(m + ra - 1) + m^ + 8mn - 9m + An^ -Qn + 2\ hm,n-i 
+ m [D(m - 2n + 1) + m^ + 12m?2 - 15m + 12n^ - 30n + 20] /i„_i,„,„i 
+ (n - 1) [/i^+2,„-2 - (£> - 3m - 4n + 4)/i^+i,„_2] • (C.l) 

D Linear Programming Implementation 

Let us write the crossing constraint Eq. (5.3) as 

= F^^{u,v) + Y^ V^,iFt:Au,v) , (D.l) 

where F^l{u,v) = v^''G/^^i{u,v) — u^"G/^^i{v,u). To rule out some spectrum of operator 
dimensions, it suffices to find a linear functional A acting on functions of {u, v) such that 

1. A(Fqq) = 1 (normalization condition) 

2. A(F^l) > for all A,/ in the spectrum (positivity constraints). 

Any such A would be inconsistent with the crossing relation Eq. (D.l) and positivity of the 
coefficients PA,h implying that the putative spectrum cannot be realized in a unitary (or 
reflection positive) CFT. 

In practice, we consider A of the form 

A:F{u,v) ^ Yl ^m,ndTd^F{a,b)\a=i,b=o (D.2) 

m+2?i<2ninax+l 
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where the variables a, b are defined in Eq. (4.15), Xm,n are real coefficients, and the range of 
m, n depends on an integer n^a.x- Since F^l{u,v) is antisymmetric under u -^ v, only odd 
a-derivatives are nonzero, and a given n^^x corresponds to (nmax + l)('^max + 2)/2 nonzero 
coefficients Xm,n- Larger n^ax gives stronger bounds, but is more computationally intensive. 
Derivatives of F^J are simply linear combinations of derivatives of the conformal blocks 
Ga,i, which we compute in Mathematica using the methods outlined in Section 4. We first 
evaluate the derivatives d^G/^j up to m = 2nmax + 1 and all the other derivatives in the 
range m + 2n < 2nmax + 1 follow via the recursion relation of Appendix C To compare with 
previous work, we have n^ax = N /2 = k — 1 where N and k are the parameters used in [17] 
and [23], respectively. 



We implement the positivity constraints above by discretizing the set of dimensions and 
restricting the spin to lie below some large finite value. The conformal blocks Ga,i{u,v) 
converge quickly for large dimensions and spins, so this is a reasonable approximation. It 
can be made arbitrarily good by using finer discretizations and a larger maximum dimension 
and spin. The plots in this paper were generated with the choices given in Table 2. 





6 


^max 


-^max 


Tl 


2 X 10-5 


3 





T2 


5 X 10"^ 


8 


6 


T3 


2 X 10"^ 


22 


20 


T4 


0.02 


100 


50 


T5 


1 


500 


100 



Table 2: In this work, we used a combination of five tables T1-T5 of conformal blocks 
(and their derivatives) with different discretizations, maximum dimensions, and maximum 
spins. For each table, dimensions were chosen from the unitarity bound Amin = I + 1 — ^6ifl 
to Amax + 2(Lmax — with step 6, and spins were restricted to < Z < Lmax- The choices 
above allow for high-resolution studies of the low-spin spectrum (T1-T3), while simultaneously 
ensuring control of intermediate dimensions and spins (T4), and also asymptotic behavior (T5). 



After restricting the dimensions and spins to lie in a finite set, our problem becomes 
a standard linear programming problem which can be solved on a computer. Solvers are 
available in a wide variety of software libraries and applications, including for example 
Mathematica. Here, we choose to use the dual simplex algorithm implementation in IBM's 
ILOG CPLEX Optimizer 
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To generate plots like those in Figure 5, we must scan over different choices of dimensions 
A„, A^, . . . , solving a linear program each time to determine the boundary between feasible 
and infeasible choices. When scanning over a single dimension, for example, this is most 
efficiently done using a binary search. One can additionally generalize binary searches to 
work in higher dimensions by recursively refining a lattice of points. These algorithms are 
readily parallelizable, and it is very convenient to take advantage of a cluster of machines 
to perform the computations. Our search logic is implemented in Scala, taking advantage 



12 



http://www-01.ibin.com/softwarc/intcgration/optiinization/cplcx-optimizcr/ 
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of its actor model for distributing parallel tasks across a network, and ILOG CPLEX's Java 
(Scala compatible) API for performing the computations. 
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